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Abstract 

The meson spectrum of QCD is studied in the framework of non- 
perturbative QCD as a function of varying quark masses m q . It is 
shown, that the total spectrum consists of two branches: 1) the stan- 
dard one, which may be called the flux-tube spectrum, depending 
approximately linearly on m q and 2) the chiral symmetry breaking 
(CSB) spectrum for pseudoscalar flavor nonsinglet (PS) mesons, with 
mass dependence yjm q . The formalism for PS mesons is derived from 
the QCD Lagrangian with m q corrections, and a unified form of the 
PS propagator was derived. It is shown, that the CSB branch of PS 
mesons joins to the flux-tube branch at around m q = 200 MeV. All 
these results are in close correspondence with recent numerical data 
on large lattices. 



1 Introduction 

The connection between the standard QCD picture of confinement, which 
may be called the flux-tube picture, and the CSB properties, which can be 
shortly called the chiral physics, is not clearly understood by physical com- 
munity, especially on the level of model building. 

The chiral Lagrangians have been introduced before the QCD era to im- 
plement CSB in [U EJ [3] and the GMOR relation |4J contains an important 
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connection between pion characteristics m n , f n and purely quark characteris- 
tics m q and (qq). From the derivation it is not clear, how this connection will 
transform with growing m q , and in general at what m q chiral mesons cease 
to be chiral and become standard flux-tube mesons, which do not contain 
chiral effects, since at large m q CSB is not spontateous any more. On the 
side of chiral Lagrangians an effective technic of Chiral Perturbation Theory 
(CPTh)[3E] was developed, which allows to calculate all corrections in terms 
of additional terms, depending on m n , f n . But one of the questions is that 
correction terms containing m q can occur already in derivation of GMOR, 
which was actually known only in the chiral limit. This will be demonstrated 
below explicitly. 

To derive GMOR and explicit expressions for m^, f n etc. we are using the 
technic of derivation of chiral degrees of freedom in QCD, developed before in 
[3 El [10]. There it is shown explicitly how pionic variables appear in QCD 
and how they are connected to the quark Green's functions. As a result 
all chiral dynamics can be expressed in terms of nonchiral (standard flux- 
tube) quark-aniquark PS Green's function G^°\p) and the so-called vertex 
(or residual) mass M(0). Both these quantities are calculated via string 
tension a and hence one can express chiral dynamics from the first principles. 

Another important point is the incorporation of the chiral branch, i.e. 
7T, K, t] and the corresponding radial excited states into the general scheme 
of mesons, or in other words, to which extent radial excitations of ir, K, e.g. 
tt(2S), K(2S) etc. are chiral objects. 

Finally, studying meson spectra with varying m q one can check the rela- 
tivistic properties of the Hamiltonian vs explicit numerical data and discover 
interesting new dependencies. Indeed, as will be shown below, the spectrum 
of vector mesons depends approximately linearly on m q with good accuracy 
both in analytic and lattice calculations, while for PS mesons linear depen- 
dence for m q > 250 MeV (flux-tube part of spectrum) goes over into the 
chiral y/m q dependence for lower m q , 




To do all analytic calculations we are using the so-called Relativistic 
String Hamiltonian (RSH) [IT] , which contains only first-principle param- 
eters (no constituent masses or subtraction constants) and was used suc- 
cessively for all hadrons: mesons [12J, baryons [13] , glueballs [H], hybrids 



Good agreement in all cases with experiment and lattice data allows to 
exploit RSH below for a precision test of flux-tube and chiral dynamics in 
QCD with varying m q . 
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The plan of the paper is as follows: in section 2 we shortly remind the 
derivation of GMOR and expressions for m n , f n in the context of QCD, and 
find additional terms in m q , which are important for growing m q . We also 
compare these results with lattice data. 

In section 3 we calculate the flux-tube spectra of vector and PS mesons 
for the whole set of m q values and compare those in detail with lattice data. 
We find here linear dependencies of flux-tube spectra. 

In section 4 we show how the chiral spectrum is incorporated in the total 
spectrum and what are chiral corrections to excited pions and kaons. Section 
5 contains summary and conclusions. 

2 Quark mass dependence of chiral dynamics 

In the standard approach one considers chiral Lagrangian and chiral per- 
turbation theory as a selfcontained approach and therefore all quark mass 
dependence (QMD) is usually derived from chiral perturbation theory. In 
what follows we shall derive chiral dynamics from QCD in a more direct way, 
as it was done in [71 El O HO], and in addition we shall calculate correction 
terms and then define QMD for all basic quantities of chiral dynamics: m n , f n 
and (ijjip). 

The starting point is the effective quark-meson Lagrangian, derived in 
[8], where both quark and chiral meson degrees of freedom (d.o.f) are kept, 

Lqml = J d A x J d 4 y{ f ip^ a (x)[{id + im f ) a p5 {4 \x - y)5 fg + 

iMU% 9 \ X) y)} ^(y) -2n f [J(x,y)]- 1 M 2 (x,y)}. (1) 

Here U^p contains Nambu-Goldstone (NG) fields <j) a , U = exp(i7 5 a t a ), 
while M(x, y) is the auxiliary scalar field. J(x, y) is expressed as an inte- 
gral of the confining field correlator D(x), its exact form will not be used 
below. 

After integrating quark d.o.f. one obtains the Effective Chiral Lagrangian 
Lecli containing effective scalar mass variable M and NG field a , with the 
partition function 

Z = I DMD<p a exp Lecl (2) 
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Figure 1: The flux-tube operator M(x,y) in the Wilson loop 



where L EC l is 

L ECL = -271// d 4 x J d 4 y(J(x,y))- l M 2 (x,y)+N c trlog{(id + im f )l + iMU]. 

(3) 

In (CQ) and (J3j) M(x, y) enters in its nonlocal form, the nonlocality being 
of the order of the vacuum correlation length A ~ 0.1 fm. In what follows 
only the local limit of M(x,y) — > M(x) will be used. In Figl. one can see, 
that M(x) can be associated with the part of the flux tube, from the quark 
position x to the center of flux tube; the same is true for antiquark. 

Expanding the last term in ([3]) in powers of <p a , one obtains the quadratic 
vertex in action 

Mr - d^k 

W^\<P) = -f] Uk)M-k)N{k)— (4) 

where 

N(k) = Ur{{AM) + J d (4 W fe2 A(0, z)M (O)A(z, 0)M(0)}. (5) 

The graphical representation of N(k) is given in Fig.2. 
Here 4> a {k) is 



f-rr fir 



jr- JL _ Si K° 

i x~ K° — / 



(6) 
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Figure 2: Quadratic terms in the pionic action (5) 



and A, A are defined as 

A= (d + m + M)" 1 , A = 0-m- M)~ l . (7) 
N(k) = hr(AMA0 - m)) = ^[tr(AM) + G (MM \k)] (8) 
and the last term in ([8]) is 

G {MM \k) = tr(AMAM) k = G {MM \k) - G (MM) (0) + G (MM) (0). (9) 

Here and in ([8]) the subscript (o, k) indicates the momentum argument. For 
A^(0) one obtains from ([8]) 

N(0) = i(tr(Am) + tr(AmAm)). (10) 

Note the last term on the r.h.s. of ffTUj) . which was not calculated before 
in [3 [U |9l [10]. Here m is the diagonal quark matrix m = diag(m u , m^, m s ), 
and we shall confine ourselves to the SU{2) case, replacing m by mu + md . 

Now trA can be expressed via the standard chiral condensate, defined in 
the Minkowskian space-time 

(^)m = i(^ + ) E = -N c trA = (uu + dd) M - (11) 
To trace down the QMD of trA it is convenient to rewrite it as 

trA = (tr(7 5 A(a;,y)75(M(0) + m)A(y, x))d 4 y = -(M(0) + m)G {0 \k = 0). 

(12) 

In (j!2p and earlier in (jSJ) we have defined the Green's functions G^ 
and G^ MM \ which differ due to presence in the latter of the operator M at 
the initial and final points x and y. In a general position M(x) contains 
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confinement interaction of the given quark with antiquark, i.e. it exemplifies 
the string between q and q, see Fig. 1. However at the initial or the end point 
of the meson propagator it can be expressed in terms of vacuum correlation 
length A [16] . A ~ 0.15 fm, namely in the Appendix 3 of [9] M(0) was 
estimated as 

M(0) = -j=<r\(l + 0(a\ 2 )). (13) 

V7T 



As shown in [8j [9j [10] and £?( MM ) can be represented as spectral sums 
over meson states 

G (0, W = -Epfr> G^)(A;) = -M 2 (0)£-|U-, (14) 
„ k z + m„ k + 

where c n = ^J^ip n (0), and ip n {x) is the n-th state meson wave function. As 
a result one obtains for (ipip) > 

m)\ = N c (M(0) +m) £ J^ffi. (15) 

n=0 m ™ 

For iV(0) from (flOl) . (|T2|) one observes a cancellation of 0(m 2 ) correction 
terms, 

mg) (Q) _ mM(0) ^ |^(0)| 2 _ mtf(Q) |(#)| 
( ) ~ 2~ (0) -^ _ ^ o ^r _ M(0)+m 2iV c (16) 

On the other hand, from the definition of N(k) in (jl]), one has a corre- 
spondence 

N(k) = « + k 2 )^- + 0(k^ (17) 
and as a consequence the (modified) GMOR relation 

mM(0) 

Ww M(0) + m IV ^ /h V ; 

which goes over into the standard GMOR formula in the chiral limit m — > 
To find fl separately one can use (JHJ), © to write 

Q( MM \k) - G (MM \0) = + 0(k 4 ) (19) 



and finally 

fl = iV c M 2 (0) £ (20) 

n=0 m n 

Eqs. (|15p .(l20 p are exploited in the Appendix to obtain numerical values 
of and fx from the calculated spectrum. 

Eqs. (ITol) . (1TF1) . (|2"U1) contain all information about QMD in the correspond- 
ing quantities: to make it explicit one needs to tell about QMD of m n and 
?/>„(0). First of all one should mention, as in [9], that all sums over n in (fT5j) . 
( !T6|) . (|20|) are formally divergent, since spectral representation fails at times 
and distances less than A, and one should use the cut-off factor in integration 
over small time region, which effectively produces the cut-off factors e~ mnX 
in summation over n in f)15p and (12"U|) . 

This fact, however, does not influence qualitatively the QMD of the cor- 
responding sums over n, which is mostly contained in the first terms of these 
sums, and appears there via m n and ip n (0, 

i*><®r = %(' + Wh ) ) (21) 

where u n = (y/p 2 + m 2 ) n . 

As will be shown in the next section by analytic calculations and com- 
parison to lattice data, the QMD of m n ,oj n is rather mild, i.e. m n grows by 
< 20% when rriQ is changing from zero to 200 MeV, and in any case this 
dependence can be computed explicitly. Neglecting it for a moment, one 
expects with this accuracy: 

a) A linear growth of with m q , i.e. for m q M(0) ~ 0.15 GeV the 
chiral condensate is twice as large as compared to zero m q value. 

b) Approximate independence of f n on m q for m q < 0.2 GeV in accor- 
dance with Eq. (1201) . 

c) Approximate scaling law mps ~ y/^h for m g < 0.2 GeV. 

More explicitly, m n and ip n (0) depend on m 2 , one can expect the following 
expansions at small m q 

m 2 PS = Xm q (l + 0(m 2 q )), (22) 
fps(m q ) = f PS (0) + O(m 2 q ), (23) 
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(i + o«)). 



(24) 



This should be compared with CPTh results [6] 




Xm q + Ym q lnm q + Zm^ 



(25) 



/ Q (l + am^ s ) + 0(m 4 PS ) 



(26) 



Note, that CPTh produces additional terms with respect to our expan- 
sions (|22p ~ (|24p . which originate from higher order terms in expansion of Lecl 
in powers of not accounted in our consideration above. 

We can now compare our expansions (12211241) with lattice data, obtained 
on large lattices 16 3 32 and 24 3 48 [17] . The detailed comparison in [17] of 
CPTh expansions f[2"5"j) . f[2"6l) with data has shown a good agreement, which 
supports also our results f[2"2"j) . (12"51) . where no serious deviations from ([25]), 
([26]) on the basis of proper quark contributions is suggested. 

The situation with the chiral condensate in quite different. Indeed, Eqs. 
([15]) and (1241) show, that (ijjip)o, entering in GMOR relation and (ipip) at 
arbitrary value of m q , differ by a linear factor K = m " ? A ^^°' > , which can be 
large, since M(0) w crA w 0(0.1 GeV) (see Appendix 1 of [9]). Additional 
corrections in m 9 , entering from the sum in f[T5]) are small for small m g < 0.2 
GeV, and the factor K has a clear significance in comparison to lattice data 
and experiment. 

In CPTh the chiral condensate is difficult to treat because of divergence 
at small m n . On lattice the condensate (qq) was studied in numerous papers, 
see e.g. [18], [19] and [20] for a review and references. Both in [18] and [19] 
the linear growth of (ipip) with m q was also observed. As one can see from 
([ED, (1201). / w ~ M(0) ~ crA, while (qq) ~ ctM(O) ~ a 2 A, where we have 
taken into account, that Vn(0) ~ uff ~ a 3 ^ 2 ,m n ~ -y/o 1 . Inserting these 
estimates into the GMOR relation (TIBl) . one obtains 



where M G ps 1/A is the gluelump mass, obtained in [16], M G ps 6-ya ~ 2 




(27) 



GeV. 
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3 The flux-tube spectrum: m ps and m v 

In the previous section the use was made of the spectrum of states {m n } and 
{ip n (0)}, which correspond to a completely different regime, described by the 
so-called Relativistic String Hamiltonian, which was derived from the path- 
integral representation of the Green's function of quark and antiquark at the 
ends of the QCD string. Therefore it exemplifies both string asymptotics at 
large L and relativistic potential dynamics at small L. For L = the RSH 
for the qq system has the form 

2 , 

H = E VP 2 + m i + V «(r) + V g (r) + A SE + H ss , (28) 

i=l 

where V a (r) = ar, V g (r) = — 4a ^ r ^ , and A S e,H ss are self-energy and 
hyperfine contributions to be defined later. 

We shall exploit the so-called einbein version of RSH, which yields results 
numerically close to (128]) . but is easier to treat. In this case one introduces 
the auxiliary variable and the total mass of the radial excited state n,n = 
0, 1, 2, ... can be written for equal quark masses m 1 = m 2 = m q , lo\ = u; 2 = to, 
as 



2 

TfL 

M n {u) = -^ + u + e n {u) + AsE{oo) + M ss (u)=Mi \u) + As E (oo) + M ss (u) 

(29) 

with 

a / 2 \ 

e n (u) = u-^a^a^r) a (0) = 2.338; r=-^f^j 



A se = -—rj{m q \), M ss (co) = ^ J fR 2 n (0)(s 1 s 2 ), (30) 



— V (m q X), M s » = ^- 2 ' 

7TU) 9iO' 

where rj is a calculable function, given in [20] note, that r](0) = 1. 

In the einbein method the equilibrium point of M^\u) is defined by the 
equation 

^° } K -o (31) 
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Table 1: Energy eigenvalues £o(m q ) and radial wave function at origin Rjs{0) 
as functions of quark mass 



m q (MeV) 





40 


100 


200 


330 


400 


500 


700 


1400 


oo ( MeV) 


352 


355 


373 


434 


507 


578 


663 


828 


1510 


T(u ),a s = 0.33 


0.3885 


0.391 


0.409 


0.447 


0.495 


0.541 


0.5925 


0.687 


1.026 


so (MeV) 


903 


900 


879 


819 


760 


711 


661 


583 


380 


R 1S (0) Gev 3 / 2 


0.312 


0.314 


0.324 


0.358 


0.397 


0.435 


0.457 


0.563 


0.911 



and Uq acquires the physical meaning of an average quark energy, ujq 1 ^ = 
m \ + P 2 )n- Note, that all correction terms, A SE and M ss , also depend on 
u, and we take them in the first approximation at the point u = Uq. 
In this way one obtains an equation for uj , 



2/3 



ul = ml + ^—a n {r) (32) 

with solutions which can be written in two forms, the first appropriate for 
large m q , 



a 




27 4 2 




(33) 



while the second is valid for small m q , when ^- < 



a 



iP I 1 
2\ — cos - arctan 



e: 

27 




(34) 



with p= q—^-- 

The resulting values of uJo(n),e r 
a = 0.18 GeV 2 and different m q . 



(ojq) are given in Table 1 for n — 0, 



For large m q , — 1 <C 1, Eq. (j32|) yields uj, 



^2 , o- 2 / 3 a n i 
" l q ^ 3 



,2/3 



One can see in Table 1, that uq change only by 3%, when m q grows from 
to 100 MeV. This fact is basic for small m q corrections to /„., (qq) from the 
sums oven in ffTS]) . f[2T?|) . 
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Table 2: Masses of vector state My and pseudoscalar state M P s together 
with spin-averaged mass M as functions of quark mass 



m q (MeV) 





40 


100 


200 


330 


400 


500 


700 


1400 


M ( MeV) 


766 


770 


874 


1001 


1224 


1357 


1535 


1890 


3154 


M v ( MeV) 


827 


830 


932 


1053 


1271 


1401 


1572 


1926 


3182 


M PS (MeV) 


458 


463 


578 


733 


984 


1135 


1349 


1709 


3011 



We now turn to corrections, A$e and M ss , which are defined as in Eq. 
( |30|) . with i?is(0) given in the Table 1, and r](m q X) in the Appendix of [20J. 
In this way one obtains for the spin-averaged mass M, at m q = 0, 

M = M cog (lS) = 766 MeV, (35) 

and for PS and V masses with cthf = 0.35 and -Ris(O) from Table 1 one would 
have 

M v = 827 MeV, M PS = 583 MeV. 

Taking now V ss interaction to the first order into account for the singlet 
wave function ijj PS (0) } one has V'ps(O) — 1-3^ps(0)' which shifts M PS down, 
M P % = 458 MeV. 

We keep this procedure for all m q and get in this way the values of Mps 
and M v , given in Table 2. 

To make comparison with lattice data we demonstrate in Table 3 numer- 
ical values of M PS , My and f P $, obtained in [17] on large lattices 16 3 32 and 
24 3 48 at (3 = 8.00 and 8.45 respectively. The corresponding values of m q were 
found using the standard procedure with chiral behavior of M PS (m q ) ~ y /m q ', 
as it is also clear from the last column of Table 3. Comparison of the values 
of My(m q ) in Table 3 and Table 2 shows an agreement within the accuracy 
of 10% 0. The same can be told about the M P s values for m q > m cr ; t ~ 200 
MeV, however for smaller m q the RSH calculation predicts M P s slowly chang- 
ing with a finite value O(400 MeV) for m q = 0. This implies that chiral 
regime is outside of RSH and around m q = 200 MeV a change of regimes 
takes place. 

^ few percent agreement occurs for a s « 0.4, but for illustrative purposes we keep in 
Tables 1 and 2a s = 0.33. 
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Table 3: Masses Mp$, My and fp$ (all in MeV), calculated in [17J, as 
functions of quark mass. 



16 3 32, 13 = 8.00 


m q (MeV) 


M PS 


My 


fps 


m £l Gev 1 / 2 


13.35 


239 


809 


94.3 


2.068 


35.2 


295 


816 


96.1 


1.57 


52.8 


353 


829 


98.1 


1.53 


70.4 


403 


841 


100.6 


1.518 


105.6 


488 


874 


103.1 


1.50 


176 


631 


944 


113.27 


1.50 


246.5 


753 


1025 


122 


1.516 


24 3 ^ 


18, (3 = 8.45 


300 


887 


1286 


153 


1.62 


211 


731 


1159 


140 


1.59 


120 


551 


1036 


125 


1.59 



One can see in Fig. 3 that our vector masses My as a function of m q with 
good accuracy lie on the sequence of two straight lines, and the RSH values 
of PS masses are on an almost parallel line for m q > m crit , while the chiral 
branch for m q < m crit follow the law Mp s ~ m q , which is in agreement with 
lattice data, shown on Fig. 3. 

In the next section we shall discuss this change of regimes from another 
side of the chiral approach and a possibility of incorporating two regimes into 
one scheme. 

It is interesting, that the entries in the rightest column of Table 3, are 
according to ( l2Tj) . the square root of the gluelump mass, y/M G , which is 
approximately equal in Table 3 to 1.5 GeV 1 / 2 , and M G ^0(2 GeV), as was 
predicted in ( j27l) . 

4 A universal spectrum in the PS channel 

The corrected GMOR relation ({TBI) gives the values of the PS masses in the 
wide interval for m q < m cr it, where the square root behavior Mps ~ 
goes over into the quasilinear RSH regime Mp$ ~ a + bm q . One might won- 
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der how the quasilinear regime, present in the spectrum of the RSH Green's 
function G^(k) = J2 n k^+m 2 * s coexisting with the chiral regime in the "to- 
tal" Green's function G(k). To this end one can compare definitions of N(k) 
in © and (O and understand, that H^N(k) = m 2 +k 2 + 0(fc 4 ) is the 

J 7T 

inverse of the generalized pion propagator, 

fl 1 
An ^ = 2N c N(k) = ml + k 2 a(k 2 ) (36) 

where 

2iV c M 2 (0) ~ cl = fflg) 



Note that a(0) = 1 due to ([H]). 

From (1361) one can find poles of A 7r (fc), which are connected to the poles 
in a(k 2 ). The latter are standard flux-tube eigenstates, the spin-averaged 
ones, since G(k) does not contain hyperfine interaction. As a consequence 



one obtains for m q = 



a(k 2 ) : poles at — k 2 = m^, n = 0,1,2, 



A w (k) : poles at k 2 = 0, ; — k 2 = m 2 n — bm 2 n (chiral), n = 1,2, ... (38) 

Note, that the m\ pole is replaced in A 7r (fc) by the pion pole (in the chiral 
limit) k 2 = 0, while m\ is shifted down by 5m\ (chiral). 

One can check, that for the initial values mo — 0.5 GeV, mi = 1.3 GeV, 
the shift of the mass eigenvalue of m n (2S) is around 0.15 GeV. This means, 
that the lowest IS eigenvalue of the RSH at m Q is replaced by m n , and the 
2S state is shifted down by ~ 10%, which is expected from the physical 
considerations, since this shift replaces the hyperfine interaction. 

In the general case of m q > the poles of A 7r (fc) are at k 2 = —M 2 ni j 

where M unif ~ ml{mq)+m 2 {mq) + U {^^))) ■ 

The situation with poles is illustrated in Fig. 4, where poles of A n (k) 
denoted as M uni f are shown together with lattice values of M Lat and the flux- 
tube values of Mps, taken from Table 2. One can see, that M uni f approaches 
Mps with growing m q and all three branches Mps, Mi at and M uni f are rather 
close to each other for m q > 150 MeV. Thus indeed, unification of chiral and 
flux-tube dynamics, automatically obtained on the lattice can be paralleled 
with our explicit mechanism of unification, given in (136)) . 
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Figure 3: Spectrum of vector and PS mesons from analytic calculations with 
the Hamiltonian fl28|) (solid and dashed lines respectively) and from lattice 
calculations [17] (dotted lines marked by the letter L). 
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Figure 4: Masses of PS mesons as functions of quark mass; M PS calculated 
in the flux-tube dynamics (solid line), M[ at , calculated on the lattice from 
[TT] . and M uni f, as poles of the generalized pion propagator, Eq. (j55|l . The 
upper solid line is for the spin-averaged masses M, from Table 2. 
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5 Summary and conclusions 



We have studied in the paper above the confining and CSB dynamics using 
varying quark masses as a tool to distinguish in the hadron spectra different 
branches, corresponding to different mechanisms. Comparison to the existing 
lattice data allows to establish a good accuracy of our relativistic formalism, 
based on RSH, and an interesting correspondence between the CSB dynamics 
and the flux-tube dynamics, and to answer the question, at what quark mass 
the Nambu-Goldstone modes transform into the standard flux-tube modes. 
In section 3 we have found, that it happens approximately at m q = m crit « 
0.2 GeV. Another interesting point is that the almost linear dependence of 

mv on m q slightly changes its slope at m' crit = \J~2o ( ^) « 0.21 GeV (cf. 
Eqs. ( 133]) and ( )34|) ). which again agrees with lattice and analytic calculations, 
as can be seen from Fig. 3. 

Maybe the most important result of the paper is the improved derivation 
of the chiral dynamics directly from the QCD Lagrangian in section 2, which 
enables one not only to obtain GMOR relation, but also find also the cor- 
responding m q dependencies. This program, started in [7J [8j |9] [10] , allows 
in principle to derive the chiral Lagrangian directly from QCD, including 
higher-order corrections, and as was already mentioned above, to connect 
chiral and flux-tube dynamics, which are usually treated separately. 

Note the important role, which is played in this program by the new entity 
- the vertex (or residual) mass M(0). Its value was calculated in terms of 
a and the correlation length A, M(0) ~ aX. Note, that in the flux-tube 
dynamics, i.e. at distances (r) « l/y/a, the mass M(0) does not appear, 
and it is needed only to calculate f % and (qq) , while in the GMOR relations 
M (0) is absent. 

Looking from the flux-tube dynamics side, it is interesting to try to under- 
stand, which kind of forces causes the breakdown of flux-tube and appearance 
of chiral dynamics, and in particular the strong reduction of the pion mass. 
To this end consider in particular, the hyperfine (hf) interaction, deduced 
in the RSH, which is singular for small values of Ui (i.e. for small average 
quark energies), since 

M ss (u) = ^^o-(s 4 s i ) 
Taking this into account and working in the einbein formalism, one can 
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find minimum of M(ui = Uj) in the PS channel, 

M(u) = l2V + ort (3g) 

CO o CO 

This minimum exists for m q > m" rit = y^o^/a ~ 0.2 GeV, but disappears 
for smaller m q , and M(u) is not bounded from below. 

Usually and in the present paper the hf interaction is treated perturba- 
tively, which in particular means, that the corresponding u in M^f is taken 
from the spin-averaged masses M^°\ 

It is not known at present how to treat M^f nonperturbatively for m q < 
m" rit . It seems very likely, that the divergence of M for small oj implies 
necessity of consideration of multiple pair creation or, in other words, the 
reconstruction of the vacuum with the appearance of quark pair condensate. 
This is still another manifestation of the appearance of the new regime - the 
chiral Nambu-Goldstone regime. 

From this point of view it is probably not surprising, that all three critical 
masses m crit ,m' crit ,m" rit approximately coincide. This fact calls for further 
studies with the aim of understanding and unifying chiral and flux-tube dy- 
namics. 

The author is grateful to A.M.Badalian for constant help and useful ad- 
vices. 



Appendix 1 

Calculation of f w and (qq) 

We illustrate in the appendix the methods of section 2 with an improved 
calculation of the chiral condensate and (qq) using Eq. f|T5l) and fl20|) re- 
spectively and the methods, given in [U [9l HO] . 

Defining these values at A -1 = 2 GeV, one should use the corresponding 
cut-off factors, and one hasl 

2 Extra factors in (jAl.ll) . (|A1.2[) . as compared to (fTij)) . ([20]) . are due to exclusion of a 
piece of the Euclidean time integration from to A, see eq. (jlll) of [9]. 
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_ M = iV c (M(0) + m q ) ^^e- m " A (Al.l) 

W / n=0 m ™ 

fl = N C M 2 (Q) £ ^V m " A (l + m n X). (A1.2) 

Here A" is taken to be 2, since higher terms are small. 

Insertion of m„, ^(0) for n = 0,1,2 from the RSH in section 3, one 
obtains 

_M^ = (217 ^V).Cj£(te) (A1.3) 
n f y ' 150 MeV v ; 

ff ^ = 96 MeV (M( ° )+mq) . (A1.4) 
^ 150 MeV V ; 

This should be compared with the values, obtained on the lattice [TO 
[T9~| [20] . which can be roughly characterized by an average value in the 
quenched case -(qq) = [(270 ± 20) MeV] 3 while from [20J for n f = 2 

(ft) _ /onnxo „„,A3 



(209 ± 8 MeV) . (Al.5) 



The results ( 1A1.1I) . ( 1A1.2I) are sensitive to the cut-off factor (vacuum 



correlation length) A, e.g. for A = 1.5 GeV 1 one obtains 

- <«>" a,V = (195 MeV) 3W)+^) (AL6) 
n f V ; 120 MeV V ; 



In the paper we adopt the results ( 1A1.3I) . (1A1.4|) and the corresponding 



estimate M(0) = 0.15 GeV as our current values. 
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